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Abstract 

When a composite of nonlinear particles suspended in a host medium is 
subjected to a sinusoidal electric field, the electrical response in the compos- 
ite will generally consist of alternating current (AC) fields at frequencies of 
higher-order harmonics. The situation becomes more interesting when the 
suspended particles are graded, with a spatial variation in the dielectric prop- 
erties. The local electric field inside the graded particles can be calculated by 
the differential effective dipole approximation, which agrees very well with a 
first-principles approach. In this work, a nonlinear differential effective dipole 
approximation and a perturbation expansion method have been employed to 
investigate the effect of gradation on the nonlinear AC responses of these com- 
posites. The results showed that the fundamental and third-harmonic AC re- 
sponses are sensitive to the dielectric-constant and/or nonlinear-susceptibility 
gradation profiles within the particles. Thus, by measuring the AC responses 
of the graded composites, it is possible to perform a real-time monitoring of 
the fabrication process of the gradation profiles within the graded particles. 
PACS number(s): 77.84.Lf, 77.22.Ej, 42.65.-k, 42.79.Ry 
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I. INTRODUCTION 



Graded materials with spatial gradients in their structure [1] are abundant in nature, 
which have received much attention as one of the advanced inhomogeneous composite ma- 
terials in diverse engineering applications [2]. These materials can be made to realize quite 
different physical properties from the homogeneous materials, and thus, to some extent, 
more useful and interesting. For graded materials, the traditional theories [3] for homo- 
geneous materials do not work any longer. Recently, we presented a first-principles ap- 
proach [4,5] and a differential effective dipole approximation (DEDA) [6,7], to investigate 
the dielectric properties of the graded materials. To our interest, the two methods have 
been demonstrated in excellent agreement between each other [4]. In the case of graded 
materials, the problem will become more complicated by the presence of nonlinearity inside 
them. Fortunately, for deriving the equivalent nonlinear susceptibility of graded particles, 
we have succeeded in putting forth a nonlinear differential effective dipole approximation 
(NDEDA) [8]. As expected, this NDEDA has also been demonstrated in excellent agreement 
with a first-principles approach [8]. 

In addition, the finite-frequency response of nonlinear composite materials have attracted 
much attention both in research and industrial applications during the last two decades [9] . 
When a composite with linear /nonlinear particles embedded in a linear /nonlinear host 
medium is subjected to a sinusoidal electric field, the electrical response in the compos- 
ite will generally consist of alternating current (AC) fields at frequencies of higher-order 
harmonics [10-14]. In fact, a convenient method of probing the nonlinear characteristics of 
the composite is to measure the harmonics of the nonlinear polarization under the applica- 
tion of a sinusoidal electric field [15]. The strength of the nonlinear polarization should be 
reflected in the magnitude of the harmonics. For the purpose of extracting such harmonics, 
the perturbation expansion [12-14] and self-consistent methods [13,16] can be used. 

In this work, based on the NDEDA, we shall investigate the effect of gradation (inho- 
mogeneity) inside the particles (inclusions) on the AC responses of the graded composite 
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by making use of a perturbation expansion method [17]. Here, the composite under con- 
sideration is composed of linear /nonlinear graded particles which are randomly embedded 
in a linear /nonlinear host medium in the dilute limit. To this end, it is shown that the 
fundamental and third-order harmonic AC responses are sensitive to the dielectric-constant 
(or nonlinear-susceptibility) gradation profile within the particle. Thus, by measuring the 
AC responses of the graded composites, it is possible to perform a real-time monitoring of 
the fabrication process of the gradation profiles within graded particles. 

This paper is organized as follows. In the following section, we shall present the formal- 
ism, which is followed by the numerical results in Section III. In Section IV, the discussion 
and conclusion will be given. 

II. FORMALISM 

Let's consider nonlinear graded spherical particles with radius a and dielectric gradation 
profile ii(r) = ei(r) + Xi( r )Ei 2 inside it, being embedded in a nonlinear host medium of 
dielectric constant e 2 = e 2 + X2E2 2 , in the presence of a uniform external electric field E 
along z— axis. Here ei(r) or e 2 (xi( r ) 01 X2) denotes the corresponding linear dielectric 
constant (nonlinear susceptibility), E x and E 2 stands for the local electric field inside the 
particles and the host medium, respectively. Note both gradation profiles ei(r) and Xi( r ) 
are radial functions where r < a. Throughout the paper, we shall consider the case of weak 
nonlinearity only, that is, Xi{ r )E\ 2 <C e\{r) and X2E2 2 <C e 2 - 

A. Comparison between a differential effective dipole approximation and a 

first-principles approach 

Recently, we put forth a DEDA (differential effective dipole approximation) [6,7] for 
calculating the equivalent dielectric constant e\{r) [8] of the spherical graded particle. This 
DEDA receives the form 
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de t (r) _ [ Cl (r)-e 1 (r)]-[e 1 (r) + 2e 1 (r)] ^ 
dr rei(r) 

Note Eq. (1) is just the Tartar formula, derived for assemblages of spheres with varying 
radial and tangential conductivity [1]. So far, the equivalent ei(r = a) for the whole graded 
particle can be calculated, at least numerically, by solving the differential equation [Eq. (1)], 
as long as ei(r) (dielectric-constant gradation profile) is given. Once e\{r = a) is determined, 
we can readily take one step forward to obtain the volume average of the linear local electric 
field inside the particles as 

< E - ( "" , > - i^fe E - < 2 > 

Hence, the DEDA [Eq. (1)] offers a convenient way to obtain the local electric field [Eq. (2)]. 
It is worth remarking that the DEDA [Eq. (1)] is valid for arbitrary gradation profiles. 

To show the correctness of Eq. (2), we shall alternatively present a first-principles ap- 
proach for calculating the local electric field inside the particle. For this purpose, let's 
take the power-law gradation profile (e(r) = A ■ (r/a) n ) as a model. For this profile, the 
potential within the graded particle can be given by solving the electrostatic equation, 
V-( ei (r)V$) = [4], 

®i(r) = -r]iE r s cos 9, r < a, (3) 

where the coefficient rji is determined by performing appropriate boundary conditions, 

3a 1_a e 2 



sA + 2e 2 ' 



and s = [V9 + 2n + n 2 — (1 + n)]/2. Based on the relation between the linear local electric 
field and the potential (Ei (lin) (r) = -V$i(r)), we have 

Ei (lin) (r) = 77i£; r s_1 {( s ~ !) cos 6» sin 6» cos 0x + (s - 1) cos 6 sin 6 sin 0y 

+ [(s-l)cos 2 fl + l]z}, (4) 

where x, y, and z are the unit vectors along x—, y—, and z— axes, respectively. So far, it is 
straightforward to obtain the volume average of the local electric field inside the particles, 



(E 1 ( Iin )) = -^E 1 ( lin )(r)dV, (5) 

where V is the volume of the spherical particles. 

In Fig. 1, we shall numerically compare Eq. (2) (local field predicted by the DEDA) with 
Eq. (5) (local field obtained from the first-principles approach). 



B. Nonlinear polarization and its higher harmonics 

1. Nonlinear differential effective dipole approximation 

In a recent work [8], we have established an NDEDA (nonlinear differential effective dipole 
approximation), by deriving a differential equation for the equivalent nonlinear susceptibility 
Xi(r), namely, 



dXi(r) - / V 

~- Xi\T) 



dr 



4dei(r)/drl _ 8y - 3 3xi(r) 
2e 2 + e!(r)J r 5r 

4 

(5 + 36:r 2 + 16:r a + 24x 4 ). (6) 



/ e 1 (r) + 2e 1 (r) , fr of; , u , , , 
, 3ei(r) 



where 



x = ei(r) - ei(r) ^ _ [e^r) - e 2 ] • [e^r) - Cl (r)] 



ei(r) + 2ei(r) ei(r)[ei(r) + 2e 2 ] 

Similarly, xx(r = a) can be obtained, at least numerically, by solving Eq. (6), once the initial 
conditions, namely, ei(r = 0) and Xiij' — 0) are given. 

In what follows, we can investigate the nonlinear AC response of the graded spherical 
particle by seeing it as a homogeneous particle having the constitutive relation between the 
displacement and the local electric field, [18] 

Di = ei(r = a)E 1 + X i(r = a)E 1 2 E 1 = l^r = a)E 1 , 

where ei(r = a) and Xi( r — a ) are determined by Eqs. (1) and (6), respectively. For the 
sake of convenience, we shall represent ei(r = a) by e±, ei(r = a) by Z\ as well as xiij' — a ) 
by Xi, if no special instructions. 
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2. Nonlinear AC responses 



If we apply a sinusoidal electric field like 

E (t) = E sm(ut), (7) 

the local electric field \J (Ei 2 ) and the induced dipole moment 

p = e e a 3 - 1 ~l 2 E (8) 
ei + 2e 2 u v ; 

will depend on time sinusoidally, too. Here the effective dielectric constant of the system 
(e e ) is given by the following dilute-limit expression, 

e e = e 2 + Se 2 f ^ jt , (9) 

6i + Ze 2 

where / is the volume fraction of the particles. By virtue of the inversion symmetry, the 
local electric field is a superposition of odd-order harmonics such that 

V^i 2 ) = E u sin(cut) + E 3ul sin(3wt) + • • • . (f 0) 

Similarly, the induced dipole moment contains harmonics as 

p = p w sin (cut) + sin(3ujt) -\ . (II) 

These harmonics coefficients can be extracted from the time dependence of the solution of 



{Ei 2 ) and p. 



3. Analytic solutions for the nonlinear AC responses 

In what follows, we will perform a perturbation expansion method to extract the third 
harmonics of the local electric field and the induced dipole moment. It is known that 
the perturbation expansion method is applicable to weak nonlinearity only, limited by the 
convergence of the series expansion. 

Let's start from the dilute-limit expression for the effective linear dielectric constant (e e ) 
of the system of interest, namely, Eq. (9) where xi — Xi — 0- 
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Next, we obtain the linear local electric fields (Ei 2 ) and {E 2 2 ), respectively, 

F 2 f)c 

(E 1 2 ) = ^^-=F(e 1 ,e 2 J,E ), (12) 
F 2 f)f 

(E 2 2 ) = -^j^=G(e 1 ,e 2 ,f,E ). (13) 

In view of the existence of nonlinearity inside the two components, we readily obtain the 
following local electric fields for the nonlinear particle and host, respectively, 

(E 1 2 ) = F^ 1 ,~e2,f,E ), (14) 
(E 2 2 ) = G(h,~t2,f,E ). (15) 

For the below series expansions, we will take e± — €± + X\E 2 ~ ei + Xi(Ei 2 ) and 
e 2 = e 2 + x 2 E 2 2 ps e 2 + X2(E 2 2 ), in Eqs. (14) and (15), where (• • •) denotes the volume 
average of • • • . Let's expand the local electric field (Ex 2 ) and (E 2 2 ) into a Taylor expansion, 
taking xi(Ei 2 ) and X2(E 2 2 ) as the perturbative quantities, 

(£i 2 > = F(e 1 ,e 2 J,E ) + ^-F(h,e 2 J, E )\ h= - ei xi{Ei 2 ) + 

J-F(e 1; e 2 , /, £ )| e ~ 2=e2 x 2 (£ 2 2 > + • • • , (16) 
(E 2 2 ) = G(e u e 2 J,E ) + ^-G(h,e 2 , f, E )\ h=Sl xi{Ei 2 ) + 

J-G(ex, e 2 , /, £ )| e ~ 2=e2 x 2 (£ 2 2 > + • • • • (17) 
Keeping the lowest orders of Xi{Ei 2 ), we can rewrite Eq. (16) as, 

(E 1 2 ) = h 1 E 2 + (h 2 + h 3 )E \ (18) 



where 



, 9e 2 2 162e 2 4 xi 
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(ex + 262) 2 ' (e! + 2e 2 )5' 

18eie 2 x 2 [(l + 3/)e! 2 + (4 - 6/)eie 2 + (4 - Qf)e 2 2 } 



(l-p)(e 1 + 2e 2 ) 5 

Because of the time-dependence of the electric field [Eq. (7)], we can take one step forward 
to obtain the local electric field in terms of the harmonics (E w and E 3ui ), 



V^i 2 ) = K sm(ut) + E 3uJ sin(3wt), 



(19) 



where 

E u = y&„ + g^^Bo 3 , (20) 

Similarly, based on Eq. (8), we obtain the induced dipole moment in terms of the har- 
monics (p w and p 3uJ ), 

P/a 3 = (p^/a 3 ) sin(wt) + (pau/a 3 ) sin(3cjt), (22) 



where 



with 



p w /a 3 = fci^b + |(fc 2 + Wo 3 , (23) 
Wa 3 = -^(*2 + W> 3 , (24) 



, ei-e 2 , 3e 2 2 /i 1 Xi[ei + 6/e 1 + (2-6/)e 2 ] 



fc 3 = 



ei + 2e 2 ' (ei + 2e 2 ) 3 

JiX 2 [(l + 3/)ei 3 - 18/ei 2 e 2 - 3(2 - 3/)e ie2 2 - 2(2 - 3/)e 2 3 ] 



(ei + 2e 2 ) 3 

In the above derivation, we have used an identity sin 3 (a>t) = (3/4) sin(cjt) — (1/4) sin(3cjt). 



III. NUMERICAL RESULTS 

For numerical calculations, we take Xi( r ) — Xi(0) + D(r/a), and other parameters: 
volume fraction p = 0.09, external field strength E = 1, linear part of host dielectric 
constant e 2 = 1. 

The validation of using the DEDA is shown in Fig. 1. In this figure, we investigate the 
linear local electric field by using a power-law gradation profile inside the particles, in an 
attempt to compare the DEDA with the first-principles approach. As expected, the excellent 
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agreement is demonstrated between the DEDA [Eq. (2)] and the first-principles approach 
[Eq. (5)]. In addition, it is worth noting that, for a linear gradation profile within the graded 
particles, the first-principles approach holds as well [4] , and the excellent agreement between 
the two methods can also be found (figure not shown here). 

Next, we discuss a power-law gradation profile [ei(r) = A ■ (r/a) n ], see Fig. 2. In this 
figure, the harmonics of local electric field and the induced dipole moment are investigated as 
a function of A for various n. In this case, increasing A (or decreasing n) leads to increasing 
ei (namely, the equivalent dielectric constant of the graded particle under consideration) and 
in turn yields a decreasing local electric field inside the particle. Thus, either an increase 
in A or a decrease in n leads to the weakening third-order harmonics (E 3oj and p 3u} ), as 
displayed in Fig. 2. 

The x— axes of Figs. 3 and 4 represent the slope (C) of a linear gradation profile. It is 
because during the fabrication of graded spherical particles by using diffusion, the dielectric 
constant at the center ei(0) may vary slightly while that at the grain boundary can change 
substantially. 

For a linear gradation profile \e\{r) — ei(0) + C ■ (r/a)], Fig. 3 shows the harmonics as 
a function of C for various ei(0). In this case, increasing C or ei(0) yield increasing e±, and 
hence one obtains the decreasing local electric field. As a result, larger C or ei(0), weaker 
the third-order harmonics (E 3uJ and p^J), see Fig. 3. 

Fig. 4 displays the effect of X2 on the harmonics, for a linear gradation profile [ei(r) = 
ei(0) + C ■ (r/a)]. Here, increasing \2 leads to increasing the local electric field inside the 
graded particle of interest. Therefore, the third-order harmonics (E 3u! and ^3^) increase for 
increasing xi- 

As mentioned above, as A and C increases, the equivalent dielectric constant of the 
particle should be increased accordingly, which in turn yields a decreasing local electric 
field, and hence, in Figs. 2~4, E u decreases for increasing A or C. On the other hand, it 
is found that, in Figs. 2~4, p^ increases for increasing A or C which is, in fact, due to the 
increasing effective dielectric constant e e [refer to Eq. (8)]. Similarly, this analysis works 
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fairly for understanding the effect of n and ei(0) on E u and p u , as displayed in Figs. 2 and 3. 
However, increasing X2 can increase not only the local electric field inside the particles, but 
also the effective dielectric constant e e , and hence we observe increasing E u and p w , as shown 
in Fig. 4. 

In addition, we also discuss the effect of nonlinear-susceptibility gradation profiles (no 
figures shown here). For linear gradation profile Xi( r ) — Xi(0) + D(r/a), as Xi(0) ( or D) 
increases, the third harmonics of both the electric field and the induced dipole moment 
increases accordingly. On the other hand, for pow-law gradation profile Xi( r ) — B(r/a) m , 
increasing B (or decreasing m) leads to increasing third harmonics. To understand such 
results, we can again resort to the above analysis on the effect of the relevant parameters 
on the local field as well as the effective dielectric constant. 

IV. DISCUSSION AND CONCLUSION 

Here some comments are in order. We investigate the nonlinear AC responses of the 
graded material where linear /nonlinear graded particles are randomly embedded in a lin- 
ear/nonlinear host medium in the dilute limit. As a matter of fact, the NDEDA (nonlinear 
differential effective dipole approximation) is valid for arbitrary gradation profiles, besides 
the power-law and linear profiles of interest. In particular, based on the first-principles ap- 
proach, the exact solution is obtainable, for not only power-law profiles (see Section II A), 
but also linear profiles (refer to Ref. [4]). 

As an extension, it is of particular interest to see what happens to the nonlinear AC 
responses of graded particles in electrorheological fluids, in which a field-induced anisotropic 
structure often occurs. For discussing this anisotropy effect, we can make use of the Maxwell- 
Garnett approximation for anisotropic structures [13]. 

To sum up, based on our recently-established NDEDA, we have investigated the nonlinear 
AC responses of a composite with linear /nonlinear graded spherical particles embedded in a 
linear /nonlinear host medium, and found the fundamental and third harmonic AC responses 
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are sensitive to the dielectric-constant (or nonlinear-susceptibility) gradation profile within 
the particles. Again, for extracting the linear local electric field, the DEDA agrees very 
well with the first-principles approach. Thus, by measuring the AC responses of the graded 
composites, it is possible to perform a real-time monitoring of the fabrication process of the 
gradation profiles within graded particles. 
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FIGURES 

FIG. 1. For power-law gradation profile ei(r) = A- (r/a) n , comparison between the approxima- 
tion result [obtained from the DEDA, Eq. (2)] and the exact solution [predicted by a first-principles 
approach, Eq. (5)], for linear electric field E±( hn ) as a function of A for various n. Parameters: 
Xi(0) = 0.1,D = 0.1, X2 = 0. 

FIG. 2. For power-law gradation profile ei(r) = A ■ (r/a) n , harmonics of the local electric field 
and induced dipole moment, as a function of A, for various n. Parameters: Xi(0) = 0.1, D = 0.1, 
X2 = 0. 

FIG. 3. For linear gradation profile ei(r) = ei(0) + C ■ (r/a), harmonics of the local electric 
field and induced dipole moment, as a function of C, for various ei(0). Parameters: xi(0) = 0.1, 
D = 0A, X 2 = 0. 



FIG. 4. Same as Fig. 3, but for various \2- Parameters: Xi(0) = 0, D = 0, ei(0) = 3. 
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